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ABSTRACT: 

We examine the status of various dark energy-models in light of the recently 
observed SN 1997ff at z ~ 1.7. The modified data still fit a pure cosmological 
constant A or a quintessence with an equation of state similar to that of A. 
The kinematical A-models, A ~ S~ 2 and A ~ H 2 , also fit the data reasonably 
well and require less dark energy density (hence more matter energy density) 
than is required by the constant A-model. However, the model A ~ S~ 2 with 
low energy density becomes unphysical as it cannot accommodate higher 
redshift objects. 

We also examine an alternative explanation of the data, viz., the absorp- 
tion by the intervening whisker-like dust and find that the quasi-steady state 
model and the FRW model f2 m o = 0.33, without any dark energy also fit the 
data reasonably well. 

We notice that the addition of SN 1997ff to the old data has worsened the 
fit to most of the models, except a closed FRW model with a constant A and 
a closed quintessence-model with = —0.82, and the models have started 
departing from each other as we go above z — 1. However, to make a clear 
discrimination possible, a few more supernovae with z > 1 are required. 

We have also calculated the age of the universe in these models and find 
that, in the models with a constant A, the expansion age is uncomfortably 
close to the age of the globular clusters. Quintessence-models show even 
lower age. The kinematical A-models are, however, interesting in this con- 
nection (especially the model A ~ H 2 ), which give remarkably large age of 
the universe. 

Subject heading: cosmology: theory, SN observations, dark energy. 
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1. INTRODUCTION 

The recent luminosity-redshift observations of type la supernovae (Riess et 
al, 1998; Perlmutter et al, 1999) reveal that the dominant constituent of 
the universe is some weird form of matter with negative pressure, commonly 
referred to as 1 dark energy 1 , which gives rise to an accelerated expansion 
and thus explains the unexpected faintness of the supernovae. This mys- 
terious cosmic matter is different from the cold, collisionless, nonbaryonic 
dark matter (CDM) (evidenced by the dynamical measurements of the to- 
tal mass density) which dominates the gravitationally clustered matter and 
clumps like the ordinary matter. The dark energy, on the contrary, does not 
clump but repels matter instead, due to its negative pressure. It is ' dark 1 be- 
cause it is not recognizable by direct observations and interacts with matter 
via gravity only. Existence of dark energy is also suggested by the anisotropy 
measurements of the microwave background radiation (MBR) which indicate 
that the universe is flat (Dodelson & Knox, 2000; Jaffe et al, 2001; Pryke 
et al, 2001; Netterfield et al, 2001). As we also have evidence that there is 
not enough matter, including dark matter, to account for the critical density 
(Turner, 2001), this provides evidence for an additional component of energy 
— the dark energy. It has been shown that the presence of dark energy in a 
CDM-model yields results that fit better with the observations of large scale 
structure (Kofman et al, 1993; Krauss & Turner, 1995; Ostriker & Stein- 
hardt, 1995; Bertschinger, 1998; Bahcall et al, 1999; Sahni & Starobinsky, 
2000). 

An obvious candidate for dark energy is a positive cosmological constant A 
which represents the energy density of the quantum vacuum (p v = A/8nG). 
However, a pure cosmological constant cannot explain the small vacuum 
energy density of 10" 120 without running into fine tuning problems. Be- 
sides this old cosmological constant problem to understand why the vacuum 
energy density is so small, now a new problem, known as the coincidence 
problem, is to understand why the vacuum energy density is comparable to 
the present mass density of the universe. That is, why we happen to live in 
this special epoch. 

Quintessence and kinematical A, which are other popular options for dark 
energy, have been proposed to solve these problems and have been extensively 
studied in the literature (for a recent review, see Sahni & Starobinsky, 2000). 
Brane-world scenarios, resulting from the hidden additional space dimen- 
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sions, have recently been proposed as an alternative to the dark energy in 
explaining the current acceleration of the universe. These theories are moti- 
vated by the fact that they are realized in string theory which is a promising 
candidate for quantum gravity. However, by considering one of such mod- 
els, it has been shown recently that the model is disfavoured by the existing 
cosmological datasets (Avelino & Martins, 2001). We hope that other brane- 
world models will be brought into the realm of cosmological testability in the 
near future. 

Supernovae la, which are thought to be thermonuclear explosions of car- 
bon oxygen white dwarfs (Trimble, 1983; Woosley & Weaver, 1986), are 
almost universally regarded as standard candles due to the following reasons. 
The variability of the luminosity of type la supernovae is small and an em- 
pirical relation exists between the decline rate of their observed luminosity 
and the absolute luminosity (de Vaucouleurs & Pence, 1976; Pskovskii 1977; 
Branch 1981, 1982; Boisseau & Wheeler 1991). That is, a brighter supernova 
fades away more slowly than the one which is less bright. Moreover, they 
have high enough absolute luminosity to make it possible to observe them at 
cosmological scales. 

As the luminosity distance c?l depends sensitively on the spatial curvature 
and the expansion dynamics of the models, the magnitude-redshift (m-z) 
relation for the high redshift supernovae has been proposed as a potential 
test for cosmological models and thus provides a useful tool to enable a 
discrimination to be made between models of various dark energy- candidates. 
Very recently, a supernova of type la (SN 1997ff) at z ~ 1.7 (the highest 
redshift SN observed so far) has been identified by the Hubble Space Telescope 
(Riess et al, 2001). It would be worthwhile to study the consequences of this 
event on various dark energy- candidates. 

In this paper, we study the consequences on various dark energy- candidates 
from the modified data, by considering the constant A FRW model and some 
variable A-models of dynamical and kinematical origin. The claim, that the 
observed dimming of the supernovae can, alternatively, be explained by the 
absorption from the whisker-like dust (Aguire, 1999; Banerjee et al, 2000), 
is also examined in the standard Friedmann model (with a vanishing A) and 
the quasi-steady state model. 

We would like to mention that there are a number of possible contam- 
inants to the measurements of SN 1997ff, for example, the host extinction 
and the magnification of the supernova by gravitational lensing due to the 
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presence of galaxies lying along our line of sight (Moertsell, Gunnarsson & 
Goobar, 2001; Lewis & Ibata, 2001). If these effects are manifested, they 
could significantly reduce the cosmological utility of this supernova, and of 
high redshift supernovae in general. However, in this analysis, we consider 
the measurements at their 'face values', which is also motivated by the tan- 
gential shear test performed and discussed by Riess et al (2001). 

Throughout our analysis, we shall be considering the Robertson- Walker 
(R-W) metric 

ds 2 = -dt 2 + S 2 (t) + r 2 (d6 2 + sin 2 Odtf)) , (1) 

[ 1 — kr 2 J 

in which case, the Einstein field equations lead to 

& k_ _ 8ttG 

g2 + ^2 ~~ Ptot, (2 J 

-g = ~y~ (Ptot + 3p to t) , (3) 

Where p t ot = Pm + Pdark energy and p t ot = Pm+Pdark energy with p m = p rcst mass + 
Pradiation Pdark matter and p m Pradiation- 



2. QUINTESSENTIAL MODEL 

It is increasingly popular to consider the possibility that the vacuum energy 
is not constant but rather evolves with the expansion of the universe. For 
example, a slowly varying scalar field (commonly known as quintessence) 
with a suitable potential V{(f>) coupled to matter through gravity, contributes 
a vacuum energy density p^ and a pressure p^ as 

p<p= 1 -4> 2 + v{4>), p = i0 2 -\/(0), (4) 

with an equation of state uj^ = p^j ' which may be a function of time 
(aty = — 1 corresponds to a pure cosmological constant). One notes that if 
the field <fi is evolving sufficiently slowly so that <fi 2 < V((f>), one can have 
uj<j, < 0. The condition (p + 3p) to t < for an accelerating expansion (see 
equation (3)) can then be satisfied. 
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As the scalar field interacts minimally with matter and radiation, these 
are conserved separately. This leads to 



p m ~ S~ 3 ^ ~ (1 + zf^\ P<p ~ S- 3 ^ 1+ ^ ~ (1 + z) 3 ( 1+w *), (5) 

where a; = and 1/3 for, respectively, rest mass (dust) and radiation, as 
usual. In deriving (5), we have assumed that is a constant which is 
a reasonable approximation as far as the equation of state changes slowly 
with time. We have also assumed that the radiation energy is negligible for 
(1 + z) « 10 3 . It is clear that p^ decays slower than p m as the potential 
energy V(<f>) starts dominating over the kinetic energy <j> 2 /2 of the scalar field 
and uj<p turns negative. However, in the earlier epoch (with p rac i ~ S* -4 ) where 
the kinetic energy is dominant, p^ decays faster (as fast as S~ & ) than p Ta £. 
Equation (5) implies that, in the present phase of evolution, 

P<p—Pm according as (l + z)—[ J , > 0, (6) 

> < \Si m o/ 

where fij's are, as usual, the energy densities of the respective components 
in units of the critical density (i.e., Jlj = 8nGpi/3H 2 ) and the subscript 
'0' denotes the value of the quantity at the present epoch. Thus in a flat 
model with Q m o = 0.34 (inferred from the SN data, as we shall see later), 
the redshift (= z cq ) at which p^ started dominating over p m , comes out as 
z eq = 0.25 if the dark energy is represented by a pure cosmological constant. 

From equations (2), (3) and (5), it is not difficult to calculate the evolution 
of the deceleration parameter q = —S/SH 2 in the present phase of evolution: 



?(*) = \ 



1 + (O^oMnoXl + 3^)(1 + z) 3 "* 



(7) 



1 + (^ /^mo)(l + Zf^ - (fi fc o/n m0 )(l + Z)- 1 

where f2 fc = k/S 2 H 2 . From this, it follows that for uo^ < —1/3 and > 0, 



q^-0 according as (1 + z) — 



(3 | lu^ I -1 



n 



(X) 



(8) 



in flat or open models. This implies that the universe must have been de- 
celerating in the past if it is accelerating at present. For example, in a flat 
universe with the dark energy in form of a constant A, the value Q m0 = 0.34 
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(which guarantees an acceleration at the present epoch), implies that the ex- 
pansion of the universe shifted from deceleration to acceleration at z — 0.57. 
This happened earlier than the epoch when the vacuum energy density took 
over to the matter energy density. 

In order to fit the models to the supernovae data, we now derive the 
magnitude (m)— redshift (z) relation. The luminosity distance of a source 
with redshift z is given by 



h{z) = ^£z[y/KH Q -£-r , (9) 




where 

£(x) = sin(x) with /C = flko when Qko > 0, 
£(x) = sinh(:r) with /C = — Qto when Qko < 0, 
£(x) = x with JC = 1 when Qko = 0. 

The apparent magnitude m of the source is, then, given by 

m( z ) =M + 51og 10 (H d L ), (10) 

where M. = M — 5 log 10 Hq + 25, M is the absolute luminosity of the source 
and o?l is measured in Mpc. It is worth noting that in deriving equations (9) 
and (10), the only assumption one needs to make, is that of the validity of the 
R-W metric given by equation (1). Hence these equations are applicable to 
all the homogeneous and isotropic models. The additional knowledge of the 
evolution of the Hubble parameter H(z) in a particular model, determines 
the m-z relation completely in the model. The expression for H(z) in this 
model follows from equations (2) and (5) as 

H(z) = H [Q m0 (l + zf + ^ (1 + zf 1+ ^ - Q fe0 (l + z?] l/2 - (11) 

We consider the data on the redshift and effective magnitude of SN 1997ff 
together with the sample of 54 supernovae from the Perlmutter et al' primary 
fit C. There are two ways Riess et al (2001) have measured SN 1997ff. One 
is to use the photometric redshift (1.7) of the supernova, in which case the 
resulting distance and redshift is a complex confidence contour (see figures 
10 and 11 in their paper). The other comes from the deployment of the 
spectroscopic redshift (1.755) of the host which has the virtue of being much 
more accurate, but is subject to being in error if the line one identifies is 
not correct. However, both methods give a very similar result and a good 
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representation of the data is z = 1.755 ± 0.05, m cS = 25.68 ± 0.34 (from a 
personal discussion with Adam G. Riess). The best-fitting parameters are 
estimated by minimizing the function 

2 _ J* mf - m(zj) 

t=l [ a mf 

where m'f refers to the effective magnitude of the ith supernovae which 
has been corrected by the lightcurve width-luminosity correction, galactic 
extinction and the K-correction from the differences of the R- and B-band 
filters and the dispersion o~ m es is the uncertainty in m| ff . 

Before fitting the data to the model, we recall that any source of energy 
and momentum should satisfy the null energy condition which, in the R-W 
metric, is equivalent to requiring p + p > 0. In the case of quintessence, this 
condition leads to p^(l +00^) > 0. The global best-fitting model, where x 2 is 
minimized by giving free rein to f2 m o, ^ao, ^4> and M., i.e. at 51 degrees of 
freedom (dof), is obtained as Q m o = 0.71, Q^o = 1.79 and uj^ = —0.82 with 
X 2 = 56.82, which is indeed a good fit (Ax 2 = x 2 / dof = L12 )- The S lobal 
best fit for a constant A {uj^ = —1) also gives a similar value, viz., x 2 — 56.89 
(at 52 dof) for the model Q m o = 0.87, Q\q = 1.51, which also shows a very 
good fit (Ax 2 = 1.09). We note that without including SN 1997ff, the global 
best-fitting quintessence model was obtained as f2 m o = 0.66, fl^o = 2.02 and 
U4> = -0.75 with x 2 = 56.81 (at 50 dof, i.e., A X 2 = 1.14). The global best- 
fitting constant A-model was: f2 m0 = 0.79, f2 A o = 1-4 with x 2 — 56.85 (at 51 
dof, i.e., Ax 2 = 1.11). 

We notice that the estimates of Q m o, mentioned above, are rather higher 
than its value estimated from various measurements, most of which seem to 
indicate a low value Q m0 « 0.33 (Turner, 2001), though, one can still argue 
that there is some more matter in an yet undetected form. Moreover, there 
are also some measurements which favour a higher value of Q m0 (Schindler, 
2001). 

If we stick to the flat models (Q m + = 1), we find that the minimum 
value of x 2 decreases as decreases. Under the constraint of the null energy 
condition, the best-fitting model is obviously that for oj^ = — 1 (i.e., for a 
constant A), which is given by fi m0 = 0.34 with x 2 — 62.01 (at 53 dof). 
This although represents a good fit (Ax 2 = 1.17), but not as good as the 
one with the non-flat models mentioned above. Nonetheless, the estimated 



(12) 
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value Q m o = 0.34 is remarkably close to the most favoured value f2 m o ~ 0.33 
from various observations, as mentioned above. The best-fitting flat model, 
without SN 1997A, is fi m0 = 0.28 with X 2 = 57.70 (at 52 dof). Thus, although 
the addition of SN 1997ff has worsened the fit to the flat model considerably, 
however, the modified data still fit either a constant A or a quintessence 
whose equation of state is similar to that of A. In the following, we shall 
limit our analysis to the flat models only, which are consistent with the 
MBR observations (which indicate that Q to t is very close to 1). 

3. PHENOMENOLOGICAL MODELS OF KINEMATIC AL A 

Models of this category, providing another popular candidate for dark energy, 
have been proposed to explain the small present value of vacuum energy 
density through a kinematical decay of A. The laws for the decay of A 
are phenomenological in nature and follow either from some dimensional 
arguments (Chen & Wu, 1990; Carvalho et al, 1992; Vishwakarma, 2000), or 
some symmetry considerations (Ozer & Taha, 1987; Vishwakarma, 2001a), 
or simply by assuming A as a function of either the cosmic time or any 
other parameter of the FRW model (for a brief summary, see, Overduin 
& Cooperstock, 1998). Here we consider the following three most popular 
models from this category: 

• case (1): A = riiS~ 2 , 

• case (2): A = n 2 H 2 , 

• case (3): A = n 3 t~ 2 , 

where the constants n^s are new cosmological parameters, replacing the orig- 
inal A, and are to be determined from the observations. The first model is 
suggested by a dimensional analysis (Chen & Wu, 1990) and also results from 
a contracted Ricci-collineation along the fluid flow vector (Vishwakarma, 
2001a) and has been the subject of most attention. The second model has 
been proposed through the similar dimensional arguments (Carvalho et al, 
1992) and also follows from other arguments. In view of the present esti- 
mates of A being of the order of Hq, this ansatz seems as a natural dynamic 
law for the decay of A. The third model has also been considered by several 
authors by imposing supplementary conditions which are equivalent to as- 
suming a power law for the scale factor (Lau, 1985; Bertolami, 1986; Berman 
& Som, 1990; Lopez & Nanopoulos 1996). However, with the flat R-W met- 
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ric, the last two models (and many other models too, for example, A ~ p, 
see, Vishwakarma, 2001b) become identical. 



Case (1): A = mS' 2 

The dynamics of the model (for k — 0) can be obtained from equations (2) 
and (3) as 

Pm = dS- 3 ^ + - ™} S \ d = constant, (13) 
47rG(l + 3w) 



^ 2 = MGd s _ (1+3w) m(l + w) 
3 " (1 + 3w) ' 

The constants n\ and C\ can be written in the following form 



(14) 



ni = 3H 2 S 2 (l-Q m0 ), d = ^^(3Q m0 - 2). (15) 

07TU 

The value of the constant C\ in the early radiation-dominated epoch can be 
calculated, in the same way, in terms of the parameters at, say, the Planck 
time or the time of the last scattering. The expression for the Hubble pa- 
rameter, in the present phase, yields from equation (14) as 

H(z) = H [(3Q m0 - 2)(1 + zf - 3(1 - fi m0 )(l + z) 2 } 1 ' 2 . (16) 

From this and equation (13), we notice that the accelerating models are 
severely constrained. As q(z) = ArcGCi/3S 3 H 2 , it is clear that for f2 m0 < 2/3 
(or Q m0 < 2f2 AO , in the non-flat models), the models have an accelerating 
expansion throughout this phase. However, the maximum redshift of an 
object in these models with a positive A, is given by z max = f2 m o/(2 — 3f2 m o) 
(in the non-flat models, z max is either 1/(2Qao — ^ m o) or Q m0 / (2Q\o — O m0 ) , 
whichever is smaller). Thus the low density models are ruled out by the 
existence of high redshift quasars. Nevertheless, the model can increase z max 
considerably by selecting f2 m o sufficiently close to 2/3 (or to 2Q\o i n the non- 
flat models), if permitted by the observations. For example, with Q m o = 0.65, 
Zmax becomes as high as 13. However, the model, then, would have to resort 
to some other explanation for the MBR, radically different from the standard 
one which assumes a cosmological origin of the MBR. 
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However no such difficulty arises for 2/3 < f2 m o < 1 (or < (H m o — 
2^ao) < 1, in the non-flat models), in which case the model is either de- 
celerating (Q m0 > 2/3) or coasting (f2 m0 = 2/3) throughout this phase of 
evolution. 

We also note that p v ^p m according as z^(2Q ia0 — l)/(2 — 3O m0 ) if fi m o > 
2/3. For Q m0 < 2/3, p v ^p m according as z^(2Q mQ — l)/(2 — 3fi m0 ). Thus, 
we find that with a positive A, the dependence of the relative dominance of 
p v and p m on z is just opposite to what we have seen in the quintessence 
model. For example, in the model Q m o = 0.65, p m started dominating over 
p v at z = 6. It should not be surprising that p v is decreasing with time faster 
than p m in this model. This is because of a continuous creation of matter 
from the decaying p v , as we shall discuss later. Equation (13), which implies 
that (p m — 2p v ) ~ (1 + z) 3 , explains why Q m and Q v are of the same order 
at the present epoch. (Note that in the coasting model, Q m = 2Q V always). 

Our fitting procedure gives the best-fitting solution as f2 m o = 0.59 with 
X 2 = 66.89 (at 53 dof). This represents a reasonably good fit (Ax 2 = 1.26), 
though not as good as we have obtained for a constant A. However, we note 
that this model can be rejected on the grounds that it cannot accommodate 
objects with redshift higher than 2.57. But the parameter space, allowed by 
the data, is wide enough to accommodate models with high enough z max . 
For example, the model f2 m0 = 0.65 gives x 2 — 68.12 (at 54 dof), which 
represents a reasonably good fit (Ax 2 = 1.26). 

Note that the best-fitting solution without including SN 1997ff is given 
by tt m0 = 0.49 with x 2 = 58.98 (at 52 dof, i.e., Ax 2 = 1.13). This implies 
that the addition of SN 1997ff has worsened the fit considerably. 

It is important to note that the quintessence fields, which acquire nega- 
tive pressure during the matter dominated phase and behave like dynamical 
A, seem formally equivalent to the kinematical models of A considered in this 
section. However, these two categories are, in general, fundamentally differ- 
ent from the theoretical point of view. In the former case, quintessence and 
matter fields are assumed to be conserved separately (through the assump- 
tion of minimal coupling of the scalar field with the matter fields) . However, 
in the latter case, the conserved quantity is [T%[ — {A(t)/87rG}g lj }, leading to 

S 

Ptot + 3(p tot + Ptot) ^ = 0, (17) 
which follows from the Bianchi identities. (Obviously, the additional require- 
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ment of the conservation of the matter field forces A to become a constant). 
This implies that there is a continuous creation of matter from the decaying 
A as is clear from the following 

o-3(l+w) , 

p m = CS- 3{1+UJ) — / k(t)S 3il+Ui) dt, C = constant, (18) 

8%G J 

which follows from (17). Thus with a positive kinematical A(t), the energy 
momentum tensor of the A-field = — (A/8rtG)g' lj (with its equation of state 
Pa = — pa = — A/8nG) acts just like the creation field of the Quasi-Steady 
State Cosmology (which will be discussed in section 6). Because of its nega- 
tive energy, the c— field also has a repulsive effect to accelerate the expansion 
as does a positive A. 

Due to this fundamental difference, the two categories of the models have, 
in general, completely different dynamics of evolution. Note that p m , in 
equations (13) and (19), no longer scales down as (1 + z) 3 ( 1+ ^) as it does 
in equation (5). In this connection, we consider an interesting example of 
the quintessential model with = —1/3 (which is the equation of state of 
the non-relativistic cosmic string). Note that the resulting model also has 
an S'" 2 -variation of the effective A, i.e., A = 8nGp ( f > ~ S~ 2 . However, the 
evolutions of p m and H obtained from equations (5) and (11), in this flat 
model, are completely different from those in (13) and (16). This is also 
obvious from figure 1 , where we have shown the expansion dynamics of these 
models, together with some other quintessential and kinematical A-models 
and compared them with the constant A-model. 
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Figure 1: The time-evolution of the scale factor S (in units of So) is 
shown as a function of the cosmic time t (in units of Hq 1 ), in some flat 
models. The plotted quasi-steady state model is given by Q\o = —0.358 
with ^ max = 5, whereas the rest of the models are specified by f2 m0 = 0.33. 



Case (2): A = n 2 H 2 

The (flat) model is described by the following equations: 

p m = c 2 S {n2 ~ z){1+w \ C 2 = constant > 0, (19) 

£ 2 = 8nGC 2 g( n2 _3)( 1+w ) +2 ^ 

(3-n 2 ) ' K } 

Note that the ansatz A = n 2 H 2 is equivalent to assuming that (= n 2 /3) 

is a constant and, hence, so is Jl m (= 1 —n 2 /3) in a flat model. We write n 2 , 
and also C 2 , in terms of the parameters at the present epoch as 

n 2 = 3(l-fi m0 ), C 2 = ^n m0 S^°. (21) 
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The expression for the Hubble parameter, in the present epoch, then yields 



H(z) = H (l + z) 3n ™» /2 . (22) 

We note that both, p m and p v scale down as (1 + z)( 3- " 2 )( 1+a; ) maintaining a 
fixed ratio p m / p v = Q m / (1 — il m ). The deceleration parameter, in the model, 
is obtained as 

t= ( 3-„Ki+*) _ lt (23) 

which is constant and implies that according as n2^(l + 3a;) /(l + a;). 
Thus two different values of the parameter 712, viz., one with n2 < 3/2 in the 
early radiation dominated era and the other with n 2 > 1 in the present era 
can make the universe shift from deceleration to acceleration. In this view, 
n 2 behaves like a tracking parameter. 

The fitting of the 54 points-data to the model gives the best-fitting solu- 
tion as Q m Q = 0.4 with x 2 — 58.35 (at 52 dof), which represents a good fit 
(A% 2 = 1.12). The addition of SN 1997ff to this dataset worsens the fit and 
the new dataset gives the best-fitting solution as f2 m o = 0.49 with x 2 = 65.27 
(at 53 dof, i.e., A% 2 = 1.23), which though represents a reasonably good fit, 
better than the one we obtained for the model A ~ S~ 2 , but not as good as 
the ones obtained for the constant A-models. 

From the best-fitting solutions of the models discussed so far (see, Table 
1), we notice that the kinematial A-models require lesser amount of dark 
energy (Qao) than is required by the quintessential or constant A-models, to 
fit the data. Said in other words, this means that with a given amount of dark 
energy density, the kinematial A-models show larger distance of an object of 
a given redshift than the distance the object has in the quintessential or 
constant A-models. This is also clear from figure 2, where we have compared 
the luminosity distance of a source (at different redshifts z) calculated in 
different models. 
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Figure 2: The luminosity distance (in units of H^ 1 ) is plotted as a 
function of the cosmological redshift in some flat models. With a given 
amount of dark energy density, the kinematial A-models show larger 
distance of an object at a given redshift than the distance the object has 
in the quintessential or constant A-models. The model A ~ S~ 2 , with 
even a smaller Q\o — 0.4, shows higher distance of a source than it has 
in the constant A- or the quintessential model with Q\o — 0.67. Note 
that all these models are specified completely by only one parameter f^o- 



It seems that all the best-fitting models discussed so far fit the data 
reasonably well. In order to be able to discriminate between the models, we 
magnify their differences by plotting the relative magnitude with respect to 
a fiducial constant A- model f2 m0 = 0.3 and Q,\0 = 0.7, as has been shown 
in figure 3. However, we realize that it is still not possible to distinguish 
the models with the available accuracy of the data. Some more high redshift 
supernovae with z > 1 are needed to make the discrimination possible, as 
the models start to depart from each other very rapidly once we go above 
z — 1, as is clear from figure 4. 



14 




Figure 3: Hubble diagram of the SN la minus a fiducial flat FRW model 
with f2 m0 = 0.3 (i.e., the relative magnitude Am(z) = m(z) — mgduciai {%))'■ 
The curves correspond to some best-fitting theoretical models. The older 
data points (z < 1) have too large error bars to make a discrimination 
possible between the models (especially between the flat models including 
the open model with zero A and whisker-dust). However, the new point 
(SN 1997ff) seems to be consistent with only the best-fitting spherical 
models. 
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Figure 4: The effective magnitude is shown as a function of the redshift 
in some best-fitting flat models and the open FRW model with zero A. 
The curves start to depart from each other very rapidly as we go above 
z = l. 

4. EFFECT OF THE INTERGALACTIC DUST 

It has been argued that the dimming of the high redshift supernovae relative 
to the low redshift ones could be an artefact of the intervening dust (Aguire, 
1999; Banerjee et al, 2000). The intergalactic dust is claimed to be domi- 
nated by needle- or whisker-type dust which is supposed to form from the 
condensation of the metallic vapours ejected from the supernovae explosions 
(Hoyle & Wickramasinghe, 1988; Narlikar et al, 1997). It is confirmed from 
the laboratory experiments that this needle-like metallic dust can grow more 
rapidly than the spherical dust (Hoyle et al, 2000). 

In a homogeneous, isotropic universe, the contribution to the effective 
magnitude arising from the absorption of light by the intervening whisker- 
like dust, is given by 
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dz' 



(24) 



H (z>y 



where k is the mass absorption coefficient, which is effectively constant over 
a wide range of wavelengths and is of the order 10 5 cm 2 g _1 (Wickramasinghe 
& Wallis, 1996), p g ~ S~ 3 is the whisker grain density with p g o in the range 
(3 — 5) x 10~ 34 g cm -3 and £(z) is the proper distance traversed by light 
through the intergalactic medium emitted at the epoch of redshift z. 
Thus the net magnitude is given by 



where the first term on the r.h.s. corresponds to the usual magnitude from 
the cosmological evolution. We consider a particular FRW model with f2 m0 = 
0.33 and A = to test equation (25). Using H Q = 70 km s~ x Mpc -1 , we find 
that the model fits the 54 point-data (i.e., without including SN 1997ff) very 
well giving x 2 = 59.78 (at 52 dof), which represents a good fit (Ax 2 = 1.15). 
However, the addition of SN 1997ff to this data worsens the fit considerably 
and the x 2 value increases to 70.38 (at 53 dof, i.e., A% 2 = 1.33), which though 
acceptable, is not as good as the ones with a non-zero A. 

Table 1. For a comparative study, the best-fitting parameters of the 
models discussed in the paper are listed, which have been estimated from 
the SN la data, with and without SN 1997ff. The parameters of all the 
models, except the QSSC, are estimated from the SN sample of Permutter 
et al' primary fit C (i.e., the excluded points are 1992bo, 1992br, 1994H, 
19970, 1996cg, 1996cn). In the case of the QSSC, the fitted data sample 
is from Banerjee et al (where the excluded points are: 1992bc, 1992br, 
1994H, 19970, 1995as, 1995K). All the other models, except the first and 
the last ones, shown in the table, are flat models. 



m net (z) =m(z) + 5m(z), 



(25) 
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0.11 


1.12 
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0.14 1.21 


0.01 


1.12 
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0.49 1.23 


0.40 


1.12 


A ~ S- 2 


0.59 1.26 


0.49 


1.13 


QSSC: z max = 5, 
Q A0 = -0.358 


1.27 




1.00 


FRW: A = 0, 
f2 m o = 0.33 


1.33 




1.15 



5. AGE CONSIDERATIONS 

The age of the universe t , in a model with the Hubble parameter H(z), is 
given by 



/•zmax H~ 1 (z) 

tn = / ±L_ ifZ dz. (26) 
Jo (1 + z) V ; 

The value of to depends on the expansion dynamics of the model considered 
and, hence, on the free parameters of the model. Additionally, the age of 
the universe depends on the present value of the Hubble parameter H . We 
consider H = 70 km s -1 Mpc -1 which comes from the present consensus 
of the Hubble Constant proponents (Freedman et al, 2001). In figure 5, we 
have plotted to as a function of fl\o in the models discussed in the preceding 
sections. We note that the age of the universe in the model with a constant 
A is uncomfortably close to the age of the globular clusters t G c = 12.5 
Gyr, even for f2 m0 as low as 0.33. The age in the quintessential models is 
even lower. Note that if the required mass density Q m o of the universe was 
smaller, one could get higher age in these models, as is clear from the figure. 
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However, this does not seem likely as the recent measurements give very 
narrow range of fl m o as Q m o = 0.330±0.035 at one sigma level (Turner, 2001). 
We also note that there are various other sources of supporting evidence that 
are consistent with this value, including studies of the evolution of cluster 
abundances with redshift, measurements of the power spectrum of large- 
scale structure, analyses of measured peculiar velocities as they relate to the 
observed matter distribution, and observations of the outflow of material 
from voids. If Q m o was, on the other hand, larger than this value, one could 
perhaps argue that there was some matter in an yet undetected form. 

In a flat model with a constant A, the fractional uncertainty in the value 
of to can be calculated, from equations (11) and (26), as 

Ato = A^o /oEWr 1 + (1 ~ Ono)^]' 372 ^ 1 - x 2 )dx A n (27) 
t H 2f 1 [n mQ x-i + (l-n m0 )x2}-^dx m0 ' 1 ' 

The recent measurement of H , by using Hubble Space Telescope, gives the 
fractional error in H as «10 percent (Freedman et al, 2001). Thus the 
expected error in t hi this model with Q m0 = 0.33 comes out as about 13 
percent (i.e., to — 13.1 Gyr, At ~ 1.7 Gyr). On the other hand, the age of 
the oldest globular cluster is estimated to be tec — 12.5 ± 1.2 Gyr (Gnedin 
et al, 2001; and the references therein). Additionally, the latest uranium 
decay estimates yield the age of the Milky Way as 12.5 ± 3 Gyr (Cayrel et 
al, 2001). Thus, if one compares the upper age limit of these objects with 
the lower limit of the age of the universe in this model, the model is in 
serious trouble. This situation in the quintessential models is even worse, as 
is clear from the figure. The situation with the best-fitting closed model with 
constant A (f2 m0 = 0.79, f2 A0 = 1.41) is not better either, which also gives 
to = 13.1 Gyr. The open model Q m0 = 0.33, A = yields t = 11.2 Gyr only. 
However, the kinematical A-models do not face this problem, especially, the 
model A ~ H 2 , which gives a remarkably large value of t = 19 Gyr in the 
best-fitting model. For a lower value of f2 m0 , to is even higher. For example, 
the value f2 m0 = 0.33 in this model gives t = 28.2 Gyr! This is not surprising 
because the expression for the age of the universe in this model yields 

which gives higher values of to for lower O m o. The model A ~ S~ 2 (with even 
a higher fl m0 = 0.65), also gives a comfortable value to = 14.7 Gyr. 
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Figure 5: The age of the universe is plotted as a function of Qao i n 
different models, by using H = 70 km s _1 Mpc -1 . The horizontal dotted 
line represents the age of the globular clusters £qc = 12.5 Gyr. The 
vertical dotted line corresponds to the mass density of the present universe 
(fi m0 = 0.33), in the flat models (Q to t = !)• In the non-flat constant 
A-model, Sl to t = 0.33 + ^aq. 



6. QUASI-STEADY STATE COSMOLOGY 

Another cosmo logical theory, where the A- energy plays an important role 
is the Quasi-Steady State Cosmology (QSSC) which is a Machian theory 
proposed by Hoyle, Burbidge and Narlikar as an alternative to the Standard 
Big Bang Cosmology (SBBC). A built-in negative A in this theory, which 
is not primarily from the observations, is necessary to make the solutions 
possible. Its role in the dynamics of the model is to energize the creation- 
field by controlling the expansion of the universe. The repulsive effects, akin 
to that from a positive A, are produced by the creation field which has a 
negative energy density. 

This cosmology does not have any cosmic epoch when the universe was 
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hot and is free from the initial singularity of the SBBC. It has cycles of 
expansion and contraction (regulated by the creation- and the negative A- 
fields respectively) of comparatively shorter period superposed on a long term 
exponential trend of expansion. Creation of matter is also periodic, being 
confined to pockets of strong gravitational fields around compact massive 
objects. The theory seems to meet all the available observational constraints 
at the present time. According to QSSC, the MBR is the relic starlight left 
by the stars of the previous cycles which has been thermalized by the whisker 
dust emitted by the supernovae. It is very interesting to note that the energy 
available from this process is just right to give a radiation background of 2. 7 
K at the present epoch (Hoyle et al, 1994). SBBC, on the other hand, does 
not predict the present temperature of the MBR. Another point in favour of 
the QSSC is that the model (which, though, has a constant A) does not face 
the cosmological constant problem mentioned in section 1 and the magnitude 
of A is broadly related with the size of the universe in a Machian way (Hoyle 
et al, 1995). In fact, the A in the QSSC does not represent the energy density 
of the quantum fields, as this model does not experience the energy scales of 
particle physics. 

In the following, we describe briefly the dynamics of the model. (For a 
detailed discussion, see Sachs et al, 1996; and the references therein). The 
field equations of QSSC, in effect, are the Einstein field equations together 
with a (negative) A and a trace-free zero rest mass scalar field c representing 
the creation of matter: 

R ij - -Rg 11 = -8ttG 

where / is a coupling constant of the c-field to gravity and q = dc/dx l . 
Assuming that the present epoch is represented by the non-creative mode of 
the model (i.e., T^L\.j = 0) and the matter is in the form of dust, equations 
(29), in the background of R-W metric, yield the following three independent 
equations: 

p m ~(l + z) 3 , c~(l + z)\ (30) 

h(z) = H [n A0 - n k0 (i + zf + n m0 (i + zf + n c0 (i + zf} 1 ' 2 , (3i) 



T ij -^9 l3 ~f[ cV + -fyg 



(29) 
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2 

[Q m0 (l + zf - 2Ov + 2ft c0 (l + z) 4 ], (32) 

where Vt c = 8nGp c /3H 2 , with p c = -3/c 2 /4 and p c = -fc 2 /A. For fc = 0, 
equation (31) reduces to 

S = S[l+r]cosilj(t)}, (33) 

where n is a parameter satisfying 0<?7<l,5isa constant given in terms 
of the constants f^o's and n and the function ip is given by 

i) 2 = -^(l + r]cosipy 2 {Q + Ar]cos^ + r] 2 (l + cos 2 ^)}. (34) 

Clearly, A has to be negative. Thus the typical energy represented by the A in 
QSSC (let us call it A-energy) has a negative energy density (= — | A | /87rG) 
and a positive pressure (=| A | /8irG). It is obvious from equation (33) that 
S never becomes zero and oscillates between 

5(1 - v ) = S min <S< S max = 5(1 + v)- (35) 

Thus at the maximum redshift z max (in the present cycle, say), we have the 
following identity from equation (31): 

^AO — ^fco(l + ^max) 2 + ^mo(l + ^max) 3 + ^c()(l + z max) 4 = 0. (36) 

As the QSSC assumes the abundance of whisker-like dust present in the 
intergalactic medium which is essential to explain the MBR, the m—z relation 
in the model is given by equation (25). Banerjee et al (2000) have already 
studied the fitting of the SN la data to this model by considering a sample 
of 54 data points which was obtained by excluding those 6 points from the 
full sample of 60 points, for which the difference between the theory and 
the observation was excessive. By assuming z max = 5, they considered a 
particular QSS model f^o = —0.358 for the actual fitting to this sample. 
The model shows an excellent fit to the sample giving x 2 = 51.97 (at 52 dof, 
i.e., Ax 2 = 1.00). The particular values of the parameters come from the 
constraints from various observations, though it is not essential to assume 
only these parameters. Indeed, the parameter space, which fits the data, is 



2q(z) 



Ho 



H(z) 
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wide enough which makes the model robust. When we add the new point 
SN 1997ff to this sample, the model still fits the resulting dataset reasonably 
well for a wide range of parameters, as is clear from figure 6, where we have 
shown a parameter-space > Q\o > —0.32 which has the A% 2 value in the 
range 1.16 — 1.26. If one sticks to the model = —0.358, then the value of 
X 2 increases to 67.41 (at 53 dof, i.e., A% 2 = 1.27), which is also acceptable. 




-0.3 -0.25 -0.2 -0.15 -0.1 -0.05 

&A0 



Figure 6: A parameter-space of interest in the QSSC: the points of the 
contour have A% 2 in the range 1.16 — 1.26. The coupling constant k and 
the present density of the whisker grains p g0 are measured, respectively, 
in units of 10 5 cm 2 g^ 1 and 10~ 34 g cm~ 3 . 



CONCLUSION 

We have studied two kinds of popular explanations of the modified SN data 
(Perlmutter st al' data together with the recently observed SN 1997ff): the 
accelerating universe-interpretation in some dark energy models of variable 
and constant A; and the absorption of light by the intergalactic dust in the 
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standard big bang model and the quasi-steady state model. The addition 
of the new point SN 1997ff to the old data has worsened the fit to most of 
the models, except a closed FRW model with a constant A, and a closed 
quintessential model with oj^ = —0.82, where it has improved the fit. In 
the flat quintessential models, the minimum value of \ 2 decreases as the 
equation of state-parameter decreases. Thus, under the constraint of the 
null energy condition, the best-fitting quintessential model has close to 
— 1. The kinenatical A-models A ~ S~ 2 and A ~ H 2 , considered in the 
paper, also fit the data reasonably well, though the fit is not as good as in 
the constant A-case. However, the model A ~ S~ 2 with low energy density 
becomes unphysical as it cannot accommodate higher redshift objects. It 
is found that the kinenatical A-models require more matter energy density 
to fit the data than is required by the constant A-model. We find that the 
expansion of the universe shifted from deceleration to acceleration at about 
z = 0.57, if the universe is flat and the dark energy is represented by a 
constant A. 

When we examine the alternative explanation of the data by considering 
the absorption of radiation by the intervening whisker-like dust, we find that 
the quasi-steady state model and the FRW model f2 m0 = 0.33, with A = 
also fit the data reasonably well, though the fit is not as good as the one to 
the constant A-model. 

While the older data was showing an almost same goodness of fit to all 
the models (see Table 1), the addition of the new point has deteriorated the 
fit to most of the models considerably and, if taken on their face values, the 
modified data show a tendency to favour a closed universe with the dark 
energy in the form of either a pure cosmological constant or a quintessence 
with uj^ = —0.82. However, the fit of the data to other models is also reason- 
ably good and some more points with z > 1 and smaller error bars would be 
needed to enable a clear discrimination between the different models. This 
might be accomplished by a dedicated space telescope SuperNova Accelera- 
tion Probe (SNAP) (Weller & Albrecht, 2001) which aims to give the accurate 
luminosity distance of ^2000 type la supernovae up to z ~ 1.7 every year. 

We have also studied the age of the universe in these models and find 
that the models with a constant A show an age which is uncomfortably close 
to the age of the globular clusters. The quintessential models yield even 
lower age. The kinematical A-models are, however, interesting in this view, 
which give large enough age, especially the model A ~ H 2 which yields a 
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remarkably large age of 28.2 Gyr for the current estimate of fi m o = 0.33. 
Note that in the QSSC, there is no beginning of the universe. 
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